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We demonstrate a controllable enhancement in the electrostrictive properties of a medium using 
dilute composite artificial materials. Analytical expressions for the composite electrostriction are 
derived and used to show that enhancement, tunability and suppression can be achieved through a 
careful choice of constituent materials. Numerical examples with Ag, AS 2 S 3 , Si and Si 02 demon¬ 
strate that even in a non-resonant regime, artificial materials can bring more than a threefold 
enhancement in the electrostriction. 
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I. INTRODUCTION 

Optoacoustic interactions have gained considerable at¬ 
tention in recent years in the context of nanophotonicsi. 
One of the strongest and most important of these is Stim¬ 
ulated Brillouin Scattering (SBS)^“— , which is a coher¬ 
ent interaction between the electromagnetic and acoustic 
fields occurring in an optical waveguide. SBS has been 
demonstrated in a number of areas within nanophotonics, 
notably in the design of nanoscale devices for Brillouin 
lasers, signal processing and microwave generation^. The 
strength of SBS is principally determined by the elec¬ 
trostriction, which is the induced strain arising from an 
electromagnetic field within the waveguiding material. 
The magnitude of the electrostrictive effect, as well as 
that of the related photo-elastic effect, has widely been 
considered a property of the material used, and as a con¬ 
sequence, the materials that have been used in SBS stud¬ 
ies have been mostly limited to those with naturally large 
electrostriction constants. 

At the same time, it is well-established in the metama¬ 
terials literature that large enhancements in the nonlin¬ 
ear properties of a medium can be achieved through the 
use of composites that have sub-wavelength structural 
features^. Metamaterials have been used to enhance non¬ 
linear scattering effects such as the Raman effect^, to 
achieve nonlinear diffractioiii, and have been used in op¬ 
tomechanical systems^ at microwave frequencies^. How¬ 
ever, nonlinear metamaterials have yet to be designed for 
the enhancement and suppression of electrostriction and 
photoelasticity, particularly in the optical range. 

In this paper we demonstrate that artificial materi¬ 
als can be designed for the tuneable enhancement or 
suppression of electrostriction. We investigate materi¬ 
als consisting of a dilute suspension of spheres embedded 
in a dielectric matrix, as presented in Fig. [TJ We con¬ 
sider both dielectric and metallic inclusions, and derive 
a mixing formula that describes the effective electrostric¬ 
tion of the composite. The electrostriction for a selection 


of practically realisable examples is then evaluated, and 
used to show that enhancement or suppression of elec¬ 
trostriction can be achieved. To our knowledge, we are 
the first to explore modifications in the optoacoustic ma¬ 
terial properties of a medium. It has been shown previ¬ 
ously that even very simple composite material designs 
can enhance the nonlinear susceptibility beyond that of 
either constituent materialai^, and therefore, we expect 
similar enhancements here with the electrostriction. 

To determine the electrostrictive properties of a com¬ 
posite material, we must first obtain electrostriction val¬ 
ues for all constituent media. Expressions for these con¬ 
stituents can differ depending on whether dispersion and 
loss are incorporated in their derivation, and less obvi¬ 
ously, on other mechanical and thermodynamic assump¬ 
tions that are imposedii^— . These considerations play 
an important part in determining regimes over which es¬ 
timates for the electrostriction are appropriate, and a 
discussion of these relevant approximations can be found 



FIG. 1. Schematic view of the metamaterial geometry inves¬ 
tigated; a primitive cubic array of spheres in a host medium. 
Inlaid: fundamental unit cell for a cubic lattice of spheres. 
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in the context of their derivations below. 

The outline of this paper is as follows. In Section [TTl 
we derive a general expression for the electrostriction, 
including the effects of dispersion, and apply this to uni¬ 
form dielectrics and metals. In Section Hill we obtain the 
electrostriction for composite materials. In Section ITVi we 
consider a series of practical examples before concluding 
remarks in Section [V] 


II. ELECTROSTRICTION FOR CONSTITUENT 
MEDIA 

In this section, we derive a general expression for the 
electrostriction of a homogeneous material. Typically, 
estimates for the electrostriction of materials are made 
under the assumption of zero loss and dispersion, zero 
shear stress, and that variations in the permittivity 
arise from changes in density alone (i.e. an isothermal 
processIn a generalisation of the standard pro¬ 
cedure, we incorporate the effects of dispersion in our 
derivation. We begin by considering the electromagnetic 
energy densityi^ 
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( 1 ) 


where Sq denotes the free-space permittivity, is the rel¬ 
ative permittivity of the material, ui is the frequency and 
E is the electric field. The change in energy density cor¬ 
responding to a small change in the density p is therefore 
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where A denotes an infinitesimal quantity. Assuming an 
isothermal process, this change in the internal energy can 
be equated to the work done W per unit volume V by 


AW = 




( 2 b) 


where P is the induced hydrostatic pressure, to obtain 


P = - 


^^o7|Ep, 


(3) 


and we define the electrostriction parameter 


d'^{£T-Oj) 

dpdoj 
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as a nondimensional measure of the induced electrostric- 
tive stress. From a microscopic perspective, the pres¬ 
sure field in ([3]) can be understood as arising from ionic 
movements in the material lattice induced by a Lorentz 
forcei^di. From we obtain expressions for the elec¬ 
trostriction of both dielectric and metallic media which 
are used in our composite model shown later. 


A. Electrostriction for dielectric media 


For a dielectric medium that is nondispersive and loss¬ 
less, the electrostriction parameter Q simplifies to the 
familiar forn:>ii 


1 = P- 


de^ 


(5a) 


It is then usual to express this in terms of other well- 
known material response tensors for practical evaluation. 
For example, for isotropic and homogeneous materials, 
(l5al) is given by 

7 = (Pii + 2pi2), (5b) 

where pij denotes the elastooptic coefficients of the 
mediuroi^ddii^. These pij coefficients are well-tabulated 
for a range of materials, and a selection of values for di¬ 
electric solids are presented in TableUfor reference. How¬ 
ever, to our knowledge no experimental data has been 
published on elastooptic coefficients for metallic mediapi^, 
and so we now consider an estimate for the 7 of metals. 


B. Electrostriction for metallic media 


For metals, we return to (|3]), and use a simplified Drude 
model for the permittivity^^ 


er = 1 - 


( 6 ) 


to obtain a form for 7 which is usefnl for practical eval¬ 
uation. Here, Wp = q‘^N/{£ome) is the square of the 
plasma frequency, q is the electric charge, rUe is the ef¬ 
fective mass of a constituent electron, N = pim is the 
number density, and m is the mass density of the metal. 
Subsequently (|3]) and (O give the estimate 


DM _ ^ 

' 9 ■ 


(7) 


We note that in the derivation of the metallic 7 above, 
we have neglected dissipation effects, which is consistent 
with the isothermal assumption made in the derivation 
of dl]). The validity of this assumption is considered in 
the results section, with a discussion of attenuation. 

Having derived evaluable expressions for dielectrics 
and metals, we now proceed to the electrostriction of 
composite materials. 


TABLE 1. Material parameters for a selection of dielectric 
materials at specified wavelengths, 7 from Eq. (I5bl) . 


Material 

A (nm) 

£r 

Pii 

P12 

7 

Ref 

Si02 

663 

2.12 

0.12 

0.27 

1.00 


AS 2 S 3 

1150 

6.06 

0.31 

0.30 

11.1 

20 

Si 

3390 

11.8 

-0.09 

0.02 

-2.77 

21 
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III. ELECTROSTRICTION FOR COMPOSITE 
MATERIALS 

In this section, we derive 7 for our composite material 
using the Maxwell-Garnett (MG) model. The effective 
permittivity given by this model is valid for a dilute ar¬ 
ray of spheres embedded in a host material, and has the 
form 2 ^ 


Since shear stress is omitted in our model, the matrix 
and inclusions do not undergo deformations, but instead 
experience a compression in order to preserve hydrostatic 
equilibrium. This is seen mathematically by stating that 
perturbations to the material pressure fields remain con¬ 
tinuous across the boundary of the sphere 


APi 


an 


AP„ 


l9n> 


( 12 ) 


_ £^m)/ 

(G + 2ein) - (G - Em)/ ' 


( 8 ) 


Here e, m denotes the relative permittivities of the con¬ 
stituent materials, and we define the filling fraction 


/ = 



(9) 


where ld,m represent corresponding volumes. The sub¬ 
script i denotes the inclusion and m denotes the matrix 
(constrained by the boundaries of the unit cell), as shown 
by the fundamental cell in Fig. [T] The only additional 
condition for ([5]) in our analysis is that the periodic array 
must be suitably subwavelength. An established rule of 
thumb is that the period of the lattice must be at least 10 
times smaller than the optical wavelength in the material. 
To demonstrate this in a practical context, a wavelength 
in the material of 1 fim and a filling fraction of / = 15% 
would correspond to an array period of 100 nm and a 
spherical radius of 33 nm. 

To begin, we consider a fully nondispersive model for 
the composite electrostriction. 


A. Nondispersive model 


Under the assumption that all constituent materials are 
nondispersive, the 7 expression o reduces to the form 
given in (I5al) . Gonsequently, from the MG model we write 


de^d^ de^ de-, de^ df 
dsm dp de-i dp df dp 


where from ([51) we have the partial derivatives 


( 10 ) 


de^ _ [(gj -b + 2/(£i - Sm)^] (1 - /) 

[(gj -b 2gin) - (gi - em)ff 
de^ _ [(gj -b 2gni) - (gj - gm)]^ / 
dSi [(£j + 2gm) - (gi - gm)/]^ ’ 

dSj. _ 3gm(gi gm)(gi T 2gjn) 

df [(gj -b 2gm) - (gi - em)ff 


( ll a) 

( ll b) 

( ll c) 


However, to evaluate the remaining three derivatives in 
m, we need to consider the mechanical response of the 
metamaterial to the induced pressure field ([3|). 


where Pi^m denotes the pressure helds and dfl is the 
boundary of the inclusion. We can then evaluate Tay¬ 
lor series for the constituent volumes Ui^m with respect 
to Pi,m to obtain 


AVi = 

AU — AP 


and express (fT^ in the form 


VA 


AW, 


an 


Kn/3n 


an 


(13a) 

(13b) 


(14) 


where we have introduced the compressibility constant 




(15) 


and K denotes the bulk modulus. Integrating both sides 
of (fT4l) we obtain the interface condition 




(16) 


for some constant A, which describes the compressive re¬ 
sponse of our composite. With this condition, and using 
the definition for the composite density 

r , n m + run, 

p = Pi/ + p4l-/) = -^F^^, (17) 


we evaluate the remaining three derivatives to obtain 


de 

m de 

m m dVi _ 7m /3m 
dp dpm dVi dp p Pc' 

dsj _ _ 7i^ 

dp dpi dVi dp p Pc' 

^ ^ ^ /(I - /) (/3m - Pi) 

dp dVi dp p Pc 


(18a) 

(18b) 

(18c) 


where analogously to (I5a1) we introduce 7 m = p-mdem/ 
dpca and 7 , = p\de\ldp\ as the electrostriction values of 
the constituent media, and Pc = /3i/-b/3m(l — /) denotes 
the volume-averaged compressibility over the unit cell. 
Gonsequently, the nondispersive electrostriction for our 
composite is given by 
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^nd 


A/ 

(Ei + 2Ei„) — ( 

^m) 

1 /?-(!-/) 

(Ei -I- 2Em)^ + 2/(Ei — Em)^ 

Pc 

(ei -|- 2Eni) — (e 

i ^m)f_ 

ft 

[(Ei -f 2Em) - (ei - £m)ff 


- A)/(l - /) 


“t” ‘^^ra) 


[(Ci + 2e^) - {£i - £m)f] 


(19) 


artificial electrostriction 


which is a weighted linear function of the constituent 
electrostriction values 71 and 7 in plus a new artificial elec¬ 
trostriction term (highlighted). The latter term can be 
understood by considering the limit 71 = 7 ni = 0 ; if the 
two materials have different compressibility values, then 
compression leads to a change in the filling fraction /, 
which, if Ei Em, alters the effective dielectric constant 
Q. Another interesting feature of (fT^ is the second- 
order pole present in all terms at 

/ = (Ci + 2Eni)/(E i — Em), (20) 

giving a theoretically infinite value for the composite elec¬ 
trostriction. However, this resonance can only be ob¬ 
tained with a change in sign for either e, or Em for dilute, 
positive /. A discussion of the asymptotic behaviour of 
(fT^ with respect to /3i,m and Ei^m is presented in Ap¬ 
pendix |A] for completeness. 

Next we consider the composite 7 expression when dis¬ 
persion is included. 


B. Dispersive corrections 


In this section, we incorporate dispersion in the deriva¬ 
tion of the composite 7 . We begin by returning to (|4]), 
which from the MG model ([ 8 ]), has the form 


i9e,. dsm de^ de^ de^ df 
_d£m dp dsi dp df dp 

d f de^ d£m de^ d£\ de^ df \ 
dp \ dsm duj dsi doj df duj j 


( 21 ) 


This composite expression is then decomposed in the 
form 7 = 7 ^^^ -|- 7 ^, where 7 ^^ and 7 ^ represent the 
nondispersive and dispersive contributions, respectively. 
The nondispersive contribution is given by the first three 
terms of m and has been evaluated in the previous 
section as (HU), where we introduce the substitution 
7i.m = therein. 

Next, we evaluate the remaining terms in (1211) . and 
note that we have 



( 22 ) 


as all mechanical parameters, such as /3 and p, are inde¬ 
pendent of the optical frequency. Accordingly, we decom¬ 
pose the dispersive term 7 ^ into a matrix and inclusion 


contribution 7 ^ = 7 ^ + 7 ° and obtain 


1m = P^- 


d£„ 


du} [(Ei -|- 2Em) ~ (£i ~ £m)/] 


^ {6EiEm/(l - /)} - ^ {6Ef/(l - /)} 
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+ ^ {ef (1 - /) - SEfEm/ - 2e^(/ + 2)- 6EiE^(l 


and 
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7^ 
Im 5 


D 9ei 
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D 


-/)} 

(23) 


(24) 


The remaining derivatives with respect to p in (1231) and 
(1241) are given in dUD where we use the substitutions 
7 i.m = 7i7 therein. 

In summary, the expression for the dispersive compos¬ 
ite electrostriction is given by 

1 = 1^^+1^+1?, (25) 

where 7 i^m = 7 i^m been substituted appropriately. 
Using the full definitions 0 for the constituents 7 i,m) 
(1^51) has an identical structure to that presented for the 
nondispersive expression in (jl9p . except now the artifi¬ 
cial electrostriction term is modified by additional terms. 
These dispersive contributions (l23l) and (|24)) feature the 
same MG resonance (uni) as before, but with a contri¬ 
bution from a third-order pole. This suggests that the 
omission of dispersion can, in certain instances, have con¬ 
siderable influence on the result for composite 7 . 


IV. NUMERICAL EXAMPLES 

In this section, we investigate the composite electrostric¬ 
tion expressions m and (1^51) for our structure, using 
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FIG. 2. Contour plots of |(a)| log^Q |Re( 7 )| from (051l . and |(b)| 
logj^g |Re(£r)| from (|8]), for an array of Ag spheres embedded 
in a Si02 matrix, against filling fraction / and wavelength A. 
A diluteness threshold of / = 15% (dashed white curve) and 
an attenuation length threshold of ql = 0 .1mm (solid white 
curve) are also shown. 


combinations of different materials. We accompany this 
investigation with an analysis of the losses for these de¬ 
signs, which is necessary for realistic applications. Ac¬ 
cordingly, we return to the MG model ([ 8 ]) and define the 
attenuation length 


ah = 


4:77 

-lm(v^) 


-1 


(26) 


which we emphasise, is completely independent from the 
electrostriction analysis. From (l26l) . a threshold of ap ^ 
0 . 1 mm is imposed as a tolerance for omitting dissipation 
effects, which is also a typical interaction length for SBS. 

We begin by investigating the composite electrostric¬ 
tion (E5|l for a cubic array of silver— spheres embedded 
in a silica^^ matrix, where we use © for 7 i. In Fig. |2(a) 


we present a contour plot of logj^p 11 ^ 6 ( 7)1 over the wave¬ 
length range 350 nm ^ A ^ 4000 nm for filling fraction 
0 0.3. 

A striking feature of this figure is the region corre¬ 
sponding to logiQ |Re( 7 )| > 0 . 8 , which contains the per¬ 
mittivity resonance (EOI)- This region simply denotes 7 
values over a cut-off threshold, which is introduced to 


ensure that features of the contour plot are not domi¬ 
nated by the singularity in (ESI) . We note that the ex¬ 
tremely strong enhancements in 7 courtesy of (I20p are 
associated with strong attenuation (E51) , and we highlight 
this by superposing a solid white curve over these con¬ 
tours, which represents an attenuation length threshold 
of ol = 0.1 mm (where to the right of this curve we have 
longer ap, and to the left, a region of shorter lengths). 

Also shown is a dashed white curve, which represents 
our diluteness threshold of / = 15%. Accordingly, in¬ 
side the region bound by these two curves (the region 
of validity (ROV)), we find a maximum composite elec¬ 
trostriction value of 7 = 3.27 at (A, /) = (1003 nm, 0.15), 
which corresponds to the intersection of the ap and / 
curves. This point gives an enhancement factor of 3.36 
relative to the electrostriction for the silica background 
at the same wavelength. It is also clear from these con¬ 
tours that the electrostriction is tuneable over a wide 
wavelength interval, but that these enhancements are ul¬ 
timately constrained by the diluteness requirement of the 
MG model. 

In Fig. |2(b)| we present a contour plot of the effec¬ 
tive permittivity over the same (A, /) range, where the 
plasmon resonance is clearly visible. For our maximum 
electrostriction value at (A, /) = (1003 nm, 0.15), we have 
a composite permittivity of = 3.4 -|- 0.003i, which is 
an enhancement factor of 1.6 relative to the background 
value of £m = 2.10 at the same wavelength. As one 
would expect, this contour plot features similar curvature 
to that of 7 , and a low degree of frequency dependence 
within the ROV (1.9 < Re(£r) < 3.4). 

We now consider silver spheres embedded in a 


we present the composite 7 for this configuration, and 
observe qualitatively similar behaviour to the previous 
example for a silica matrix in Fig. |3(a)[ The primary 
difference here is the much more restrictive ol thresh¬ 
old, which now extends to much longer wavelengths. 

If one searches inside the ROV constrained by our ol 
and / bounds, we discover a maximum electrostriction 
value of 7 = 27.4 at (A, /) = (2064 nm, 0.15), correspond¬ 
ing to the intersection of the ol and / curves as before, 
wit h an enhancement factor of 2.63 (c.f., 7 m = 10.44). 
Fig |3(b)| reveals that this coordinate point has an effec¬ 
tive permittivity value of £r = 9.3 -I- O.Oli. This corre¬ 
sponds to a similar permittivity enhancement factor as 
the previous example (c.f., £m = 5.89). A slightly higher 
level of frequency dependence is observed in the ROV 
also (5.8 < Re(£r) < 9.3). 

For these examples, we find that the composite elec¬ 
trostriction expression (l25l) gives a 10 — 20% increase 
in the maximum electrostriction value compared to the 
nondispersive expression dl). This suggests that the 
omission of dispersion can give rise to a small but non- 
negligible correction to the composite electrostriction. 
Furthermore, a similar investigation with An spheres em¬ 
bedded in these matrix materials reveals a comparable 
level of enhancement to Ag. 


chalcogenide^S matrix (amorphous AS 2 S 3 ). In Fig. 3(a' 
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FIG. 3. Contour plots of |(a)| log^Q |Re( 7 )| from (051l . and |(b)| 
logj^g |Re(£r)| from (| 8 ]), for an array of Ag spheres embedded 
in a AS 2 S 3 matrix, against filling fraction / and wavelength 
A. A diluteness threshold of / = 15% (dashed white curve) 
and an attenuation length threshold of ctl = 0 . 1 mm (solid 
white curve) are also shown. 


FIG. 4. Contour plots of |(a)| log^g |Re( 7 )| from (IT^ . and |(b)| 
logjo |Re(er)| from (| 8 ]), for an array of Si 02 spheres embedded 
in a Si matrix, against filling fraction / and wavelength A. A 
diluteness threshold of / = 15% (dashed white curve) and 
an attenuation length threshold of ol = 0 . 1 mm (solid white 
curve) are also shown. 


In Fig. |4(a)| we consider log^Q |Re( 7 )| from (fTOl) for sil¬ 
ica spheres embedded in a silicon^l matrix. This figure 
exhibits strong frequency dependence for A < 1000 nm 
(courtesy of a material resonance for Si at A ~ 370 nm) 
and a near-horizontal arc of zero electrostriction which 
spans the entire ROV. That is, this metamaterial design 
can completely suppress electrostriction over an excep¬ 
tionally wide frequency range. For this particular com¬ 
posite the attenuation length threshold is reached at ap¬ 
proximately A = 1000 nm. In Fig. |4(b)| we present a 
contour plot of log^Q |Re(er.)| for completeness, which ex¬ 
hibits reassuringly minimal frequency dependence over 
the ROV. 

In Fig. |5(a)| we show a cross section of the composite 
Re( 7 ) from Fig. |4(a)| at A = 1550 nm. This gives con¬ 
firmation that complete suppression of electrostriction 
is achieved at / ~ 10 %, and shows that we have sign¬ 
changing electrostriction from this metamaterial design. 
We note that the composite Re( 7 ) (blue curve) exceeds 
that of the constituent electrostriction values (dashed 
curves) at a filling fraction of / = 16.6%, which is rem¬ 
iniscent of earlier work which showed the nonlinear pa¬ 


rameters of composite materials can exceed the values of 
the constituents^^, but we note that care must be taken 
as the dilute lattice assumption breaks down in this re¬ 
gion of enhancement here. An investigation using AS 2 S 3 
spheres in a Si matrix demonstrated an identical result 
to that shown in Figures 4(a) 4(b) and |5(a)[ but at much 
lower filling fractions. 

To emphasise the result presented in Fig. 5(a)[ we 
show the composite 7 curve for an array of chalcogenide 
(AS 2 S 3 ) sp heres embedded in a silica matrix at A = 1550 
nm in Fig. |5(b)[ This shows a simple linear enhancement 
from the background electrostriction, to a maximum re¬ 
alisable value of 7 = 1.918 at the threshold of / = 15% 
(i.e. an enhancement factor of approximately 2 ). 


V. CONCLUDING REMARKS 


We have presented an analytical representation for the 
electrostriction of a composite material by incorporat¬ 
ing the simplest and analytically most transparent model 
from effective index theory, the MaxwelRGarnett model. 
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(a) 



(b) 


serial description of all stress fields, and is the next step 
in the development of the theory. Also, we note that 
our estimate for the electrostriction of metals is a low- 
order approximation, which requires experimental data 
for validation. 

As a final comment, we emphasise that other ho¬ 
mogenisation procedures'^ can be used to determine the 
electrostriction for a periodic composite, which should re¬ 
move several constraints of the present MG model, and 
open the way to investigations of exciting metamaterial 
designs. 

Note added in proof. Recently, we became aware of a 
preprint^ on a similar topic. 
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Appendix A: Asymptotic analysis of the 
nondispersive composite electrostriction 

In this appendix we examine several asymptotic limits 
for the composite electrostriction expression (HU). Eirst, 
assuming e„i 3> ei, we obtain 


FIG. 5. Composite electrostriction Re( 7 ) from (1191) at A = 
1550 nm for an array of (a) Si 02 spheres embedded in a Si 
matrix, and (b) AS 2 S 3 spheres embedded in a Si 02 matrix. 
The electrostriction for the inclusion 71 (dashed black line) 
and matrix 7 m (dashed red line) materials are also given. 


to the problem of electrostriction. 

We show that expressions for the electrostriction of 
a composite material feature artificial electrostriction 
terms, which contribute to the enhancement or sup¬ 
pression of this material property, as observed for a se¬ 
lection of composites here. The presence of this term 
points towards the possibility that large enhancements 
in 7 , beyond both material values, could be achieved 
for more sophisticated metamaterial designs. We also 
show that sign-switching electrostriction is achievable, 
and that resonant enhancements in the electrostriction 
of metal-dielectric composites are unrealistic, as they are 
associated with strong attenuation. Incorporating dis¬ 
persive effects in the model is shown to give a small but 
non-negligible correction to estimates for the composite 
electrostriction. 

It is important to emphasise that this work is a first 
step in the study of the electrostriction of composites, 
and so other considerations such as thermally-induced 
electrostriction and scattering losses, are not addressed 
here. These effects are more prominent for high-intensity 
wave problems, where more elaborate models are re¬ 
quired to accurately evaluate the material response. In¬ 
cluding the effect of shear stresses will require a fully ten- 


W 6(A - /3m)/(l - /)em 

/3c(/ + 2)2^' /3c(/ + 2)^”^ /3c(/ + 2)2 

(Al) 


where we have a persistent, but simplified, contribution 
from all terms in (TT^ . By contrast, for £i ^ Em we have 

9/A /3m(l + 2/) - A) 

^ ^ (1 - /)2a V £?) (1 - /)/3c (1 - me 

(A2) 

The limit An ^ A gives the asymptotic form 


(ei + 2£i„)^ + 2/(ei — Em)^ 

[(Ei -I- 2Em) - (ei - £ni)ff 


+ 


[(Ei -f 2Em) - (C - £m)ff 


, (A3) 


and A ^ /3m leads to 


7 ^ 


(Ei 2Em) (^i ^m) 

(^i T 2Eni) (e, £^m)/_ 

_ 3(1 - /)Em(Ei - 
[(ei -|- 2Em) ~ 


7i 


^m)(^i T 2Eiji) 

{Si - £m)f]^ 


(A4) 


These difference in sign in the artificial electrostriction 
contributions above suggest that the relative magnitudes 
of /3 are relevant in establishing whether enhanced or 
suppressed electrostriction is observed. 
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